Estimating Sensor Motion from Wide-Field Optical Flow on a Log-Dipolar Sensor

Log-polar image architectures, motivated by the structure of the human visual field, have long been investigated in computer vision for use in estimating motion parameters from an optical flow vector field. Practical problems with this approach have been: (i) dependence on assumed alignment of the visual and motion axes; (ii) sensitivity to occlusion form moving and stationary objects in the central visual field, where much of the numerical sensitivity is concentrated; and (iii) inaccuracy of the log-polar architecture (which is an approximation to the central 20°) for wide-field biological vision. In the present paper, we show that an algorithm based on generalization of the log-polar architecture; termed the log-dipolar sensor, provides a large improvement in performance relative to the usual log-polar sampling. Specifically, our algorithm: (i) is tolerant of large misalignmnet of the optical and motion axes; (ii) is insensitive to significant occlusion by objects of unknown motion; and (iii) represents a more correct analogy to the wide-field structure of human vision. Using the Helmholtz-Hodge decomposition to estimate the optical flow vector field on a log-dipolar sensor, we demonstrate these advantages, using synthetic optical flow maps as well as natural image sequences

Exact Geosedics and Shortest Paths on Polyhedral Surface

We present two algorithms for computing distances along a non-convex polyhedral surface. The first algorithm computes exact minimal-geodesic distances and the second algorithm combines these distances to compute exact shortest-path distances along the surface. Both algorithms have been extended to compute the exact minimalgeodesic paths and shortest paths. These algorithms have been implemented and validated on surfaces for which the correct solutions are known, in order to verify the accuracy and to measure the run-time performance, which is cubic or less for each algorithm. The exact-distance computations carried out by these algorithms are feasible for large-scale surfaces containing tens of thousands of vertices, and are a necessary component of near-isometric surface flattening methods that accurately transform curved manifolds into flat representations.National Institute for Biomedical Imaging and Bioengineering (R01 EB001550

Multi-Area Visuotopic Map Complexes in Macaque Striate and Extra-striate Cortex

We propose that a simple, closed-form mathematical expression--the Wedge-Dipole mapping--provides a concise approximation to the full-field, two-dimensional topographic structure of macaque V1, V2, and V3. A single map function, which we term a map complex, acts as a simultaneous descriptor of all three areas. Quantitative estimation of the Wedge-Dipole parameters is provided via 2DG data of central-field V1 topography and a publicly available data set of full-field macaque V1 and V2 topography. Good quantitative agreement is obtained between the data and the model presented here. The increasing importance of fMRI-based brain imaging motivates the development of more sophisticated two-dimensional models of cortical visuotopy, in contrast to the one-dimensional approximations that have been in common use. One reason is that topography has traditionally supplied an important aspect of "ground truth", or validation, for brain imaging, suggesting that further development of high-resolution fMRI will be facilitated by this data analysis. In addition, several important insights into the nature of cortical topography follows from this work. The presence of anisotropy in cortical magnification factor is shown to follow mathematically from the shared boundary conditions at the V1-V2 and V2-V3 borders, and therefore may not causally follow from the existence of columnar systems in these areas, as is widely assumed. An application of the Wedge-Dipole model to localizing aspects of visual processing to specific cortical areas--extending previous work in correlating V1 cortical magnification factor to retinal anatomy or visual psychophysics data--is briefly discussed.National Institute of Health/National Institute of Biomedical Imaging and Bioengineering (R01 EB001550

Sparsity-Promoting Calibration for GRAPPA Accelerated Parallel MRI Reconstruction

The amount of calibration data needed to produce images of adequate quality can prevent auto-calibrating parallel imaging reconstruction methods like generalized autocalibrating partially parallel acquisitions (GRAPPA) from achieving a high total acceleration factor. To improve the quality of calibration when the number of auto-calibration signal (ACS) lines is restricted, we propose a sparsity-promoting regularized calibration method that finds a GRAPPA kernel consistent with the ACS fit equations that yields jointly sparse reconstructed coil channel images. Several experiments evaluate the performance of the proposed method relative to unregularized and existing regularized calibration methods for both low-quality and underdetermined fits from the ACS lines. These experiments demonstrate that the proposed method, like other regularization methods, is capable of mitigating noise amplification, and in addition, the proposed method is particularly effective at minimizing coherent aliasing artifacts caused by poor kernel calibration in real data. Using the proposed method, we can increase the total achievable acceleration while reducing degradation of the reconstructed image better than existing regularized calibration methods.National Science Foundation (U.S.) (CAREER Grant 0643836)National Institutes of Health (U.S.) (Grant NIH R01 EB007942)National Institutes of Health (U.S.) (Grant NIH R01 EB006847)National Institutes of Health (U.S.) (Grant NIH P41 RR014075)National Institutes of Health (U.S.) (Grant NIH K01 EB011498)National Institutes of Health (U.S.) (Grant NIH F32 EB015914)National Science Foundation (U.S.). Graduate Research Fellowship Progra

Accelerated parallel magnetic resonance imaging reconstruction using joint estimation with a sparse signal model

Accelerating magnetic resonance imaging (MRI) by reducing the number of acquired k-space scan lines benefits conventional MRI significantly by decreasing the time subjects remain in the magnet. In this paper, we formulate a novel method for Joint estimation from Undersampled LinEs in Parallel MRI (JULEP) that simultaneously calibrates the GeneRalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) reconstruction kernel and reconstructs the full multi-channel k-space. We employ a joint sparsity signal model for the channel images in conjunction with observation models for both the acquired data and GRAPPA reconstructed k-space. We demonstrate using real MRI data that JULEP outperforms conventional GRAPPA reconstruction at high levels of undersampling, increasing the peak-signal-to-noise ratio by up to 10 dB.National Science Foundation (U.S.) (CAREER Grant 0643836)National Center for Research Resources (U.S.) (P41 RR014075)National Institutes of Health (U.S.) (NIH R01 EB007942)National Institutes of Health (U.S.) (NIH R01 EB006847)Siemens CorporationNational Science Foundation (U.S.). Graduate Research Fellowship Progra

Intersubject Regularity in the Intrinsic Shape of Human V1

Previous studies have reported considerable intersubject variability in the three-dimensional geometry of the human primary visual cortex (V1). Here we demonstrate that much of this variability is due to extrinsic geometric features of the cortical folds, and that the intrinsic shape of V1 is similar across individuals. V1 was imaged in ten ex vivo human hemispheres using high-resolution (200 μm) structural magnetic resonance imaging at high field strength (7 T). Manual tracings of the stria of Gennari were used to construct a surface representation, which was computationally flattened into the plane with minimal metric distortion. The instrinsic shape of V1 was determined from the boundary of the planar representation of the stria. An ellipse provided a simple parametric shape model that was a good approximation to the boundary of flattened V1. The aspect ration of the best-fitting ellipse was found to be consistent across subject, with a mean of 1.85 and standard deviation of 0.12. Optimal rigid alignment of size-normalized V1 produced greater overlap than that achieved by previous studies using different registration methods. A shape analysis of published macaque data indicated that the intrinsic shape of macaque V1 is also stereotyped, and similar to the human V1 shape. Previoud measurements of the functional boundary of V1 in human and macaque are in close agreement with these results

Ultra-high spatial resolution BOLD fMRI in humans using combined segmented-accelerated VFA-FLEET with a recursive RF pulse design

Purpose To alleviate the spatial encoding limitations of single-shot EPI by developing multi-shot segmented EPI for ultra-high-resolution fMRI with reduced ghosting artifacts from subject motion and respiration. Methods Segmented EPI can reduce readout duration and reduce acceleration factors, however, the time elapsed between segment acquisitions (on the order of seconds) can result in intermittent ghosting, limiting its use for fMRI. Here, "FLEET" segment ordering--where segments are looped over before slices--was combined with a variable flip angle progression (VFA-FLEET) to improve inter-segment fidelity and maximize signal for fMRI. Scaling a sinc pulse's flip angle for each segment (VFA-FLEET-Sinc) produced inconsistent slice profiles and ghosting, therefore, a recursive Shinnar-Le Roux (SLR) RF pulse design was developed (VFA-FLEET-SLR) to generate unique pulses for every segment that together produce consistent slice profiles and signals. Results The temporal stability of VFA-FLEET-SLR was compared against conventional-segmented EPI and VFA-FLEET-Sinc at 3 T and 7 T. VFA-FLEET-SLR showed reductions in both intermittent and stable ghosting compared to conventional-segmented and VFA-FLEET-Sinc, resulting in improved image quality with a minor trade-off in temporal SNR. Combining VFA-FLEET-SLR with acceleration, we achieved a 0.6-mm isotropic acquisition at 7 T--without zoomed imaging or partial Fourier--demonstrating reliable detection of BOLD responses to a visual stimulus. To counteract the increased repetition time from segmentation, simultaneous multi-slice VFA-FLEET-SLR was demonstrated using RF-encoded controlled aliasing. Conclusions VFA-FLEET with a recursive RF pulse design supports acquisitions with low levels of artifact and spatial blur, enabling fMRI at previously inaccessible spatial resolutions with a "full-brain" field of view.Comment: 51 pages (including supplement), 8 main figures, 6 supporting figures. For supporting videos (8), please visit https://github.com/aveberman/vfa-fleet. Note: this work has been accepted for publication at Magnetic Resonance in Medicin

Activation of the CREB/c-Fos pathway during long-term synaptic plasticity in the cerebellum granular layer

The induction of long-term potentiation and depression (LTP and LTD) is thought to trigger gene expression and protein synthesis, leading to consolidation of synaptic and neuronal changes. However, while LTP and LTD have been proposed to play important roles for sensori-motor learning in the cerebellum granular layer, their association with these mechanisms remained unclear. Here, we have investigated phosphorylation of the cAMP-responsive element binding protein (CREB) and activation of the immediate early gene c-Fos pathway following the induction of synaptic plasticity by thetaburst stimulation (TBS) in acute cerebellar slices. LTP and LTD were localized using voltage-sensitive dye imaging (VSDi). At two time points following TBS (15 min and 120 min), corresponding to the early and late phases of plasticity, slices were fixed and processed to evaluate CREB phosphorylation (P-CREB) and c-FOS protein levels, as well as Creb and c-Fos mRNA expression. High levels of P-CREB and Creb/c-Fos were detected before those of c-FOS, as expected if CREB phosphorylation triggered gene expression followed by protein synthesis. No differences between control slices and slices stimulated with TBS were observed in the presence of an N-methyl-Daspartate receptor (NMDAR) antagonist. Interestingly, activation of the CREB/c-Fos system showed a relevant degree of colocalization with long-term synaptic plasticity. These results show that NMDAR-dependent plasticity at the cerebellum input stage bears about transcriptional and post-transcriptional processes potentially contributing to cerebellar learning and memory consolidation